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DEFENCE  RESEARCH  ESTABLISHMENT  SUFFIELD 
RALSTON  ALBERTA 


SUFFIELD  TECHNICAL  NOTE  NO.  435 


The  density,  vapor  pressure,  surface  tension  and  viscosity 
of  thickened  Dowanol,  DPM,  were  measured  experimentally  according  to 
standard  techniques.  Each  of  the  techniques  is  explained  in  detail  and 
a brief  summary  of  the  theory  is  also  presented. 


A new  parameter  called  the  ^dissolved  density*  was  defined 
which  promises  to  allow  differentiation  between  thickening  polymers. 

The  addition  of  an  acrylate  polymer  had  no  effect  on  the  surface  tension 
or  vapor  pressure  of  the  solvent.  The  effect  of  thickener  concentration 
on  viscosity  was  determined  over  a range  of  shear  rates  where  the  thic- 
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SUFFIELD  TECHNICAL  NOTE  NO.  435 


PHYSICAL  PROPERTIES  OF  A THICKENED 


SOLVENT  (U) 


C.L.  Chenier  and  G.A.  Hill 


INTRODUCTION 


There  are  several  important  physical  properties  that  can  be 
measured  to  determine  the  characteristics  of  materials.  Primary  among 
these  are  tin  chemical  purity  and  toxicity.  The  vapor  hazard  and  flow 
properties  are  also  frequently  important.  In  order  to  quantify  these 
last  two  parameters,  such  things  as  density,  vapor  pressure,  surface 
tension  and  viscosity  are  commonly  measured. 


In  chemical  defence  training  exercises,  solvents  are  commonly 
employed  which  have  low  toxicity  properties  and  can  be  detected  on 
three-way  detector  paper+  or  by  standard  analytical  techniques.  One  such 


solvent  is  Dowanol,  DPMTT  (di propylene  glycol  methyl  ether).  The  manu- 


t Federal  Stock  Number  FSN  6665-21-0309 

ft  Dowanol,  DPM  is  a trade  name  of  the  Dow  Chemical  Company. 
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facturer  of  this  solvent  provides  exhaustive  data  on  most  of  its  chemical, 
physical  and  toxicological  properties. 

When  small  amounts  of  Rohm  and  Haas  K-125  (an  acrylate  thickener) 
are  added,  the  viscosity  of  Dowanol  is  observed  to  increase.  This  report 
provides  quantitative  information  on  this  viscosity  change  as  well  as 
the  effect  on  density,  surface  tension  and  vapor  pressure. 

EXPERIMENTAL  PROCEDURES 

Most  of  the  techniques  used  in  these  experiments  are  standard 
and  available  in  many  laboratory  manuals  or  chemistry  textbooks  (see  Refs. 
1,  2 and  3).  They  are  included  in  Appendix  I of  this  report  for  con- 
venience and  easy  reference. 

Several  standard,  Newtonian  fluids  were  chosen  for  calibrating 
the  equipment  used  in  the  above  procedures.  These  fluids  are  commercially 
available  and  are  listed  in  the  Tables  as  Standards  A,  B,  C,  etc.  The 
fluid  to  be  tested  was  Dowanol,  DPM,  a commercial  solvent  available  only  in 
45  gallon  drums.  Gas  chromatography  showed  our  sample  to  be  88%  pure,  the 
other  12%  being  basically  propylene  glycol  methyl  ether  and  tripropylene 
glycol  methyl  ether.  The  thickener  was  Rohm  and  Haas  K-125,  a commercial 
powder  consisting  of  80%  methyl  methacrylate  and  20%  butyl  and  ethyl  acry- 
lates. 


THEORY 

The  theories  behind  the  techniques  used  in  these  experiments 
are  well  described  in  most  books  on  physical  chemistry  (e.g.  Ref.  3).  A 
good  book  on  the  field  of  rheology  is  Reference  4.  In  the  following 
paragraphs,  formulas  will  be  presented  with  little  or  no  derivations. 

a.  Density  (see  Fig.  1) 


P 


s 


M p 
s w 


(Eq.  1) 
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where:  ps  = density  of  unknown  species  (g/cm3) 

M$  = mass  of  species  in  pycnometer  (g) 

M = mass  of  water  in  pycnometer  (g) 
w 

p = density  of  water  (g/cm3) 

W 

Equation  1 is  used  to  calculate  the  density  of  the  polymer 
solutions.  If  a high  density  substance  is  mixed  with  a low  density  sub- 
stance, the  density  of  the  resulting  material  will  lie  between  the  densities 
of  the  starting  materials  (assuming  no  chemical  reactions).  Therefore,  if 
the  density  of  the  solution  increases  with  the  addition  of  the  polymer,  the 
polymer  can  be  assumed  to  have  a higher  density  than  the  pure  liquid. 
Mathematically  this  is  expressed  by  the  following  equation. 

mpps  = Pp^mp  + mL  " VLps^  (E(l- 

where:  mp  = mass  of  polymer  (g) 

p$  = density  of  solution  (g/cm3  ...  from  Eq.  1) 

Pp  = density  of  polymer  (g/cm3) 
mL  = mass  of  pure  liquid  (g) 

VL  = volume  of  pure  liquid  (cm3) 

Equation  2 can  be  used  to  evaluate  the  dissolved  density  of  the  polymer 
by  plotting  (mpPs)  vs  (mp  + mL  - Vlp$). 

b.  Surface  Tension  (see  Fig.  2) 


Ys  = 


rp$  gh 


(Eq.  3) 


where: 


•y  = surface  tension  of  species  (dynes/cm) 


r = radius  of  capillary  (cm) 
g = gravitational  acceleration  (981  cm/sec2) 
h = length  of  rise  of  meniscus  (cm) 

c.  Vapor  Pressure  (see  Fig.  3) 

This  parameter  was  measured  directly  with  mercury  manometers 
and  a precision  mercury,  glass  thermometer.  With  pressure  and  temperature 
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data,  it  is  possible  to  calculate  the  heat  of  vaporization  from  the 
Clausius-Clapeyron  equation: 

log  P = - AHvap/(2. 303RT)  + Constant  (Eq.  4) 

where:  P = measured  pressure  (dynes/cm2) 

T = measured  temperature  (°K) 

R = gas  constant  (8.314  x 107  ergs/g  mole/°K) 

AHvap  = ^eat  vaP°rizati°n  (ergs/g  mole) 

d.  Ostwald  Viscosity  (see  Fig.  4) 

n = Kp$t  (Eq.  5) 

where:  n = shearing  viscosity  (g/cm/sec) 

K = instrument  constant  (cm2/sec2) 
t = measured  time  interval  (>  200  sec) 

The  shear  rate  cannot  be  controlled  with  these  viscometers  but  depends  on 
the  capillary  diameter  and  the  viscosity  of  the  fluid  being  used.  Also, 
because  the  pressure  head  varies  during  the  experiment,  the  shear  rate  is 
not  constant  throughout  the  experiment.  An  average,  Newtonian  shear  rate 
can  be  calculated  for  each  experiment  from  the  following  formula: 

y = 4V/Ur3t)  (Eq.  6) 

where:  V = volume  of  bulb  (3.2  cm3) 

r = radius  of  capillary  (cm) 

y = shear  rate  (sec-1) 

e.  Brookfield  Viscosity  (see  Figs.  5 and  6) 

Using  this  instrument,  the  shearing  viscosity  can  be  determined 
as  a function  of  shear  rate.  The  shear  rate  for  a Newtonian  fluid 
is  (Ref.  4): 

Y = 2n/(l  - s2)  (Eq.  7) 

where:  s = R^/Rq  (dimensionless) 

Ri  = radius  of  spindle  (cm) 

Rq  = radius  of  housing  cylinder  (cm) 
fi  = rotation  rate  (sec-1) 
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For  a non-Newtonian  fluid  this  expression  becomes  slightly  more  complex 
(see  Appendix  II): 


y = 2 n l s2Pm  m 

P=o 


(Eq.  8) 


where: 


m 


Ri 


slope  of  log  n vs  log  tRi- 

evaluated  at  t = s2PtRi-  (dimensionless) 

shear  stress  on  spindle  (dynes/cm2) 


xRi  = T/(2*R.2Le) 


The  shear  stress  is  calculated  from: 

(Eq.  9) 

where:  T = measured  torque  on  inner  cylinder  (dyne-cm) 

Le  = equivalent  length  of  spindle  (cm) 

Because  the  inner  spindle  has  surface  area  at  both  top  and  bottom, 
viscous  drag  occurs  not  only  between  the  two  cylindrical  faces  but  also 
at  the  top  and  bottom  of  the  spindle.  This  effect  is  compensated  for 
by  using  an  equivalent  length  for  the  spindle  instead  of  its  actual 
length.  This  equivalent  length  is  calculated  by  using  a standard,  New- 
tonian fluid  of  known  viscosity  and  using  the  formula: 

(Eq.  10) 

Once  y is  determined  from  equation  6 or  7 and  tRi-  is  determined  from 
equation  8,  the  shearing  viscosity  is  calculated  from: 


Le  = T/(27rR.2nY) 


n = tr./y 


(Eq.  11) 


Unfortunately,  for  non-Newtonian  liquids  n is  a function  of  y,  so  when 
xR.j  is  plotted  against  y a straight  line  is  not  obtained.  The  type  of 
curve  which  results,  however,  helps  to  characterize  the  fluid  as  being  a 
Bingham  plastic,  shear-thinning  or  shear-thickening. 

Generally,  non-Newtonian  behavior  is  modelled  by  the  simple 
power-law  equation: 


Ri 


V 


(Eq.  12) 


where:  nQ  = shear  viscosity  when  y = 1 sec"1 

n = power- law  coefficient 
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This  study  did  not  investigate  time  dependent  behavior  or 
elastic  properties  of  the  fluids. 

RESULTS 

The  data  obtained  from  these  experiments  are  listed  in  Tables  1 
to  10  and  are  depicted  in  Figures  7 to  21  for  ease  of  visualization. 

DISCUSSION  OF  THE  RESULTS 

a.  Density 

As  expected,  the  densities  of  the  polymer  solutions  were  greater 
than  the  density  of  the  pure  solvent.  The  density  increased  linearly  with 
the  increase  in  concentration  as  shown  by  Fig.  7.  An  equation  relating 
density  to  polymer  concentration  for  thickened  Dowanol  is: 

P$  = 0.959  + 0.000199c 

where:  p$  = solution  density  (g/cm3) 

c = polymer  concentration  (g/s,) 

Fig.  8 shows  that  the  relation  represented  by  Equation  2 is 
indeed  linear.  The  slope  of  the  line  is  equal  to  the  dissolved  density 
of  the  polymer.  For  K-125  in  Dowanol  the  dissolved  density  is  1.209  g/cm3. 
This  number  could  be  used  for  characterizing  unknown  polymers  or  deter- 
mining the  similarity  between  batches  of  the  "same"  polymer.  This  number 
might  not  be  identical  for  different  solvents  because  the  dissolved 
density  will  be  a function  of  the  molecular  association  between  the  solvent 
and  the  polymer. 


b.  Surface  Tension 

As  seen  in  Table  2,  the  surface  tension  of  Dowanol  remained 
constant  at  29.3  ±0.2  dynes/cm  and  was  not  affected  by  the  addition  of 
K-125. 


c.  Vapor  Pressure 

Table  4 and  Figure  9 present  data  for  two  runs  with  unthickened 
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Dowanol  and  one  run  with  thickened  Dowanol  (DO-4)  on  the  Ramsey-Young 
apparatus.  In  Figure  9 it  can  be  clearly  seen  that  within  the  experi- 
mental accuracy  no  difference  could  be  measured  between  the  vapor  pres- 
sures of  thickened  and  unthickened  Dowanol.  A least  squares  line  through 
the  data  yielded  the  following  equation  for  vapor  pressure: 

log  P = - 2457/T  + 7.942 

where:  P = kilopascals 

T = °K 

The  accuracies  on  the  constants  at  the  90%  confidence  limits  are: 

- 2457  ± 275 
and  7.942  ± 0.106 

Experimental  inaccuracies  appeared  to  be  due  to  the  impurities  in  the 
original  Dowanol  samples. 

From  the  value  of  the  slope.  Equation  4 can  be  used  to  calcu- 
late the  corresponding  heat  of  vaporization  (AH  ).  For  Dowanol,  DPM, 

AHvap  is  4.7  x 1011  ergs/g  mole. 

d.  Equivalent  Lengths 

The  five  Newtonian  fluids  listed  as  standards  in  Table  1 were 
used  to  evaluate  the  equivalent  lengths  of  the  spindles  used  in  the  Brook- 
field Small  Sample  Adaptor.  These  values  are  listed  in  Table  5 after 
calculation  by  Equation  10.  It  should  be  noted  that  even  though  spindle  3 
had  an  actual  length  longer  than  spindle  1,  its  equivalent  length  was 
significantly  smaller.  The  actual  length  to  gap  ratio  for  spindle  3 was 
five  times  larger  than  that  of  spindle  1,  thus  the  end  effects  were  more 
significant  for  spindle  1 than  for  spindle  3.  This  accounts  for  the 
greater  equivalent  length  of  spindle  1.  Both  spindles  1 and  2 had  equiva- 
lent lengths  significantly  greater  than  their  actual  lengths,  showing  that 
end  effects  were  contributing  to  the  viscous  drag  on  these  spindles. 

e.  Viscosities 

(1)  Qstwald  - The  Ostwald  experimental  viscosities  calculated 
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by  Equation  5 are  listed  in  Table  2 along  with  the  corresponding  average 
shear  rates  (calculated  by  Equation  6).  The  viscosity  is  seen  to  increase 
with  polymer  concentration.  However,  because  the  shear  rate  differs  for 
each  viscosity  determination,  it  would  not  be  proper  to  compare  these 
numbers  quantitatively.  For  non-Newtonian  fluids  the  viscosity  changes 
with  the  shear  rate  and  this  factor  could  influence  any  comparison  made 
from  the  Ostwald  viscometers.  It  is  necessary  to  compare  the  viscosities 
at  identical  shear  rates  and  the  Ostwald  viscometers  are  not  capable  of 
doing  this. 

(2)  Brookfield  - This  instrument  gives  experimental  data  re- 
lating revolutions  per  minute  to  torque.  Using  Equation  9 of  the  Theory, 
this  data  is  converted  to  revolutions  per  minute  versus  shear  stress  as 
shown  in  Table  6.  One  needs  to  convert  the  RPM  data  to  shear  rate  data. 
There  are  two  alternative  methods  of  doing  this.  Equation  7 can  be  used 
to  calculate  the  shear  rate  directly.  However,  this  assumes  the  liquid 
is  Newtonian.  The  other  method  is  to  use  Equation  8 which  takes  into 
account  the  non-Newtonian  characteristics  of  the  fluid.  Because  extremely 
viscous,  non-Newtonian  fluids  were  being  used,  both  techniques  were  tried. 
To  use  Equation  8,  one  must  first  evaluate  'm',  the  slope  of  log  TRi  vs 
log  a.  Least  squares  linear  regression  lines  were  calculated  for  the  data 
of  Table  6 and  yielded  the  values  of  m shown  in  Table  7.  The  fact  that  m 
differs  significantly  from  1.0  for  some  of  the  more  viscous  liquids  demon- 
strates the  non-Newtonian  behavior  of  these  liquids.  Using  these  values 
of  m,  the  true  shear  rates  for  each  fluid  were  evaluated  and  are  listed  in 
Table  8 along  with  the  approximate  Newtonian  shear  rates.  Although  values 
of  m as  high  as  1.36  were  found,  the  Newtonian  shear  rates  varied  at  most 
only  2%  from  the  true  shear  rates.  This  small  deviation  cannot  justify 
the  extra  calculations.  Figures  10  to  18  demonstrate  that  there  is 
essentally  no  difference  between  the  non-Newtonian  and  Newtonian  plots 
of  shear  stress  vs  shear  rate.  Thus,  it  is  concluded  that  the  Newtonian 
shear  rate  formula  is  accurate  enough  .or  fluids  whose  'm'  coefficient 
is  less  than  1.4. 
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From  the  linear  regression  lines  shown  in  Figures  10  to  18, 
the  viscosities  at  y = 1.0  sec'1  and  the  power  law  coefficients  are  easily 
evaluated  from  the  intercept  and  slope  data.  These  values  are  listed  in 
Table  9.  Thus,  for  example,  the  power  law  equation  relating  shear  stress 
to  shear  rate  for  DO-7  is: 

t = (6.08  ± 0.42)  y (°-882  1 °-021) 


The  error  limits  for  these  constants  are  based  on  the  90%  confidence 
limits  described  by  Neville  and  Kennedy  (Ref.  5).  It  should  be  mentioned 
that  these  equations  are  valid  only  over  the  limited  range  of  shear  rates 
shown  in  Table  8 for  which  the  data  were  taken. 


A quantitative  comparison  of  viscosity  can  be  made  for  each 
thickened  fluid.  Figure  19  shows  the  viscosity  ratio  (viscosity  of  solu- 
tion divided  by  viscosity  of  solvent  ...  both  at  y = 1 sec-1)  as  a 
function  of  the  polymer  concentration.  The  curve  is  similar  to  those 
shown  by  Middleman  (Ref.  4)  and  at  high  viscosity  ratios  (>  20)  it  can  be 
seen  that: 


,r  - c4-9  ...  (V20)  .... 


where  nr  = viscosity  ratio  (dimensionless) 


c = polymer  concentration  (g/2.) 

This  relation  agrees  closely  with  that  observed  by  many  others  (Ref.  4). 


Finally,  in  Table  10  and  Figures  20  and  21,  the  effect  of  shear 
rate  on  viscosity  is  demonstrated  for  all  the  thickened  fluids.  The  curves 
shown  on  Figures  20  and  21  represent  the  theoretical  relationship  given  by: 


n = 


where  n is  the  power  law  coefficient  romputed  and  listed  in  Table  9.  We 
see  that  the  data  and  the  curvec  on  Figures  20  and  21  show  a shear  thinning 
effect  for  the  high  concentration  solutions  and  essentially  a Newtonian 
relationship  for  the  low  concentration  solutions.  Figures  20  and  21  also 
demonstrate  that  the  power  law  equation  is  inaccurate  at  very  low  shear 
rates  where  it  predicts  the  viscosity  will  approach  infinity.  The  Ellis 
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model  could  be  used  to  overcome  this  problem  but  was  avoided  here  because 
it  necessitates  the  determination  of  an  extra  constant  (Ref.  4). 

CONCLUSIONS 

1.  The  density  of  thickened  Dowanol  is  a linear  function  of  the 
concentration  of  polymer  K125  up  to  50  g/t.  The  dissolved  density  of 
polymer  K125  in  Dowanol  is  1.209  g/cm3. 

2.  The  surface  tension  is  unaffected  by  polymer  addition. 

3.  The  vapor  pressure  is  unaffected  by  polymer  addition. 

The  heat  of  vaporization  of  Dowanol  is  4.7  x 10n  ergs/g  mole. 

4.  Although  Dowanol  is  a Newtonian  liquid,  the  addition  of  polymer 
K125  produced  non-Newtonian,  shear  thinning  solutions. 

5.  The  relationship  between  log  viscosity  ratio  vs  log  polymer 
concentration  is  linear  at  high  viscosity  ratios  with  a slope  of  5. 

6.  The  viscosity  of  all  solutions  could  be  modelled  by  the  power 
law  equation  over  shear  rates  of  0.5  sec-1  to  500  sec'1. 
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TABLE  1:  FLUIDS  L 

SED 

Fluid  Description 

Report  Acronym 

Brookfield  Standard  100 

ST-A 

Brookfield  Standard  1000 

ST-B 

ASTM  Oil  Standard,  S20 

ST-C 

ASTM  Oil  Standard,  S200 

ST-D 

ASTM  Oil  Standard,  S600 

ST-E 

0.0  g K-125/litre  Dowanol 

DO-1 

9.92  g K-125/litre  Dowanol 

DO-2 

20.00  g K-l 25/litre  Dowanol 

DO-3 

30.05  g K-l 25/ litre  Dowanol 

DO-4 

30.13  g K-l 25/ litre  Dowanol 

DO- 5 

36.61  g K-125/litre  Dowanol 

DO-6 

39.94  g K-l 25/litre  Dowanol 

DO- 7 

f . 49.97  g K-125/litre  Dowanol 

DO-8 

50.20  g K-l 25/litre  Dowanol 

DO-9 

UNCLASSIFIED 

4 

UNCLASSIFIED 


TABLE  2:  DENSITIES,  SURFACE  TENSIONS  AND 
OSTWALD  VISCOSITIES  AT  25°C 


Liquid 

Acronymn 

Density 

(g/cm3) 

Surface  Tension 
(dynes/cm) 

Ostwald  Viscosity 
(poise) 

ST-A 

0.9626 

20.55 

0.949 

ST-B 

0.9699 

21.54 

9.92 

ST-C 

0.8669 

30.31 

0.3313 

ST-D 

0.8758 

31.87 

4.073 

ST-E 

0.8834 

32.90 

12.46 

DO-1 

0.9505 

29.36 

0.03585 

(576)+ 

DO- 2 

0.9524 

29.22 

0.1396 

(168) 

DO- 3 

0.9537 

29.18 

0.4683 

(64) 

DO-4 

0.9562 

29.36 

1.612 

(24) 

DO- 5 

0.9570 

29.20 

1.556 

(29) 

DO-6 

0.9573 

29.78 

3.475 

(16) 

DO- 7 

0.9577 

29.32 

5.144 

(13) 

DO-8 

0.9612 

29.45 

18.32 

(3) 

DO- 9 

0.9599 

i 

28.93 

18.41 

(3) 

+ Numbers  in  brackets  represent  average  shear  rates  (sec-1) 
calculated  by  Equation  6. 
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TABLE  3:  DISSOLVED  DENSITY  OF  POLYMER  DATA 


Dowanol 

Sample 

mp  * ”L  - Vs+ 

(g) 

Vs+ 

(g2/cm3) 

DO-1 

0.00 

0.00 

DO-2 

8.02 

9.45 

DO- 3 

16.80 

19.07 

DO-4 

24.35 

28.73 

DO- 5 

23.63 

28.83 

DO- 6 

29.80 

35.05 

DO- 7 

32.74 

38.25 

DO- 8 

38.27 

48.03 

DO-9 

40.80 

48.19 

+ 

for  an  explanation  of  variables  see  Equation  2 ... 
data  plotted  as  X vs  Y in  Figure  8 


f . 
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TABLE  4:  VAPOR  PRESSURE  DATA 


DO-1 

DO-1 

DO-4 

Temp  (°K) 

Pres  (kPa) 

Temp  (°K) 

Pres  (kPa) 

Temp  (°K) 

Pres  (kPa) 

411.95 

91.86 

413.25 

84.53 

413.95 

93.46 

403.15 

87.86 

408.25 

78.79 

412.15 

89.73 

401.95 

82.66 

404.35 

73.19 

409.75 

84.26 

399.65 

77.06 

401.45 

68.13 

407.55 

78.66 

398.75 

71.59 

401.35 

63.06 

405.15 

73.33 

395.05 

65.86 

396.75 

57.46 

403.65 

67.46 

392.55 

60.26 

394.25 

52.13 

400.35 

61.86 

391.75 

54.53 

390.25 

46.80 

397.95 

56.26 

388.35 

49.06 

389.95 

41.20 

395.35 

50.66 

386.85 

43.06 

388.95 

35.86 

391 . 55 

44.80 

383.55 

38.00 

387.55 

30.80 

390.15 

39.46 

379.75 

32.40 

380.95 

25.06 

385.75 

33.86 

376.45 

26.80 

376.15 

20.00 

379.65 

28.13 

368.65 

21.06 

369.15 

14.53 

374.45 

22.53 

361.15 

15.60 

354.75 

9.47 

366.45 

17.07 

349.85 

10.00 

353.15 

6.93 

354.55 

11.47 

335.15 

4.40 

341.35 

5.87 
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TABLE  5:  EQUIVALENT  LENGTHS  OF 
BROOKFIELD  SPINDLES 


Spindle 

Acronym 

Actual  Length 
(cm) 

Equivalent  Length 
(cm) 

SP-1 

0.8382 

3.150 

3.52  ± 0.03 

SP-2 

0.8382 

6.299 

6.58  ± 0.05 

SP-3 

0.9322 

3.162 

3.19  ± 0.08 

Equivalent  Lengths  Determined  With  Three  Standard  Liquids 
Radius  of  Cylinder  Housing  - 0.9525  cm 


r. 
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TABLE  6:  SHEAR  STRESS  DATA  FOR  DOWANOL 


tr1  (dynes/cm2) 


V 


Spindle 

(1 

Sample 

(acronymn) 

(sec'1) 

00-1 

DO- 2 

DO- 3 

DO-4 

00-5 

DO-6 

DO-7 

DO-8 

DO-9 

0.0524 

. 

- 

- 

- 

. 

. 

- 

- 

- 

0.1047 

- 

- 

- 

- 

- 

- 

4.632 

- 

18.06 

0.2618 

- 

- 

- 

- 

- 

8.801 

12.51 

44.46 

43.54 

SP-1 

0.5236 

- 

- 

- 

7.874 

7.411 

17.14 

24.55 

85.22 

82.45 

1.047 

- 

1.390 

4.632 

15.75 

14.82 

33.35 

47.71 

148.2 

148.2 

2.094 

- 

2.779 

9.264 

30.11 

29.18 

63.00 

88.93 

255.7 

258.0 

5.236 

1.390 

7.411 

22.23 

69.48 

69.02 

140.4 

193.2 

503.9 

500.2 

10.47 

2.779 

12.97 

42.61 

130.2 

127.8 

248.7 

334.9 

807.8 

807.8 

0.0524 

- 

. 

- 

. 

- 

- 

8.906 

0.1047 

- 

- 

- 

- 

- 

- 

5.195 

- 

17.81 

0.2618 

- 

- 

- 

- 

- 

8.906 

12.37 

- 

43.40 

SP-2 

0.5236 

- 

- 

- 

7.422 

- 

17.32 

24.49 

- 

81.64 

1.047 

- 

1.484 

4.701 

14.84 

- 

33.15 

46.76 

- 

159.1 

2.094 

- 

2.969 

9.154 

29.19 

- 

62.59 

87.83 

- 

251.4 

5.236 

1.732 

6.680 

22.27 

68.53 

- 

139.3 

190.2 

- 

492.9 

10.47 

3.464 

13.36 

42.30 

126.4 

- 

246.6 

322.6 

- 

801.6 

0.0524 

- 

- 

. 

- 

9.906 

12.79 

46.48 

47.46 

0.1047 

- 

- 

8.667 

7.885 

18.57 

26.41 

86.32 

87.50 

0.2618 

- 

- 

- 

20.64 

19.92 

44.16 

62.32 

179.3 

186.5 

SP-3 

0.5236 

- 

- 

11.97 

40.45 

39.01 

82.55 

114.3 

305.4 

316.6 

1.047 

- 

6.604 

23.52 

75.12 

73.87 

149.0 

201.4 

507.9 

538.2 

2.094 

3.302 

14.03 

45.81 

137.0 

135.3 

260.0 

347.5 

820.0 

861.8 

5.236 

9.080 

34.26 

105.2 

296.3 

293.0 

521.7 

683.5 

1504 

1516 

10.47 

17.75 

66.45 

195.2 

531.6 

521.6 

924.5 

1179 

2350 

2371 
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TABLE  7:  DETERMINATION  OF  m 


Spindle 

Acronymn 

Fluid 

Acronymn 

Slope1- 

X 

Deviation1'1' 

dx 

Coefficient 

m = 1/x 

Accuracy1'1' 
dm  = dx/x2 

1 4 2 

DO-1 

1.000 

±0.657 

1.000 

±0.657 

1 4 2 

DO-2 

0.964 

±0.036 

1.037 

±0.039 

1 4 2 

DO- 3 

0.959 

±0.009 

1.043 

±0.010 

1 4 2 

DO-4 

0.940 

±0.018 

1.064 

±0.020 

1 

DO- 5 

0.951 

±0.028 

1.052 

±0.031 

1 4 2 

DO-6 

0.904 

±0.014 

1.106 

±0.017 

1 4 2 

DO- 7 

0.914 

±0.023 

1.094 

±0.028 

1 

DO-8 

0.781 

±0.042 

1.280 

±0.069 

1 4 2 

DO-9 

0.840 

±0.031 

1.190 

±0.044 

3 

DO-1 

1.048 

±0.242 

0.954 

±0.220 

3 

DO-2 

0.999 

±0.057 

1.001 

±0.057 

3 

DO- 3 

0.931 

±0.033 

1.074 

±0.038 

3 

DO-4 

0.892 

±0.024 

1.121 

±0.030 

3 

DO-5 

0.923 

±0.042 

1.083 

±0.049 

3 

DO-6 

0.855 

±0.028 

1.170 

±0.038 

3 

DO-7 

0.844 

±0.039 

1.185 

±0.055 

3 

DO-8 

0.736 

±0.031 

1.359 

±0.057 

3 

DO-9 

0.735 

±0.038 

1.361 

±0.070 

+ Slope  of  log  rRi  vs  log  ft 
++  90%  confidence  limit 
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TABLE  9:  SHEAR  VISCOSITY  AT  y - 1.0  AND 
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TABLE  10:  SHEAR  RATE  EFFECT  ON  VISCOSITY 


Newtonian  Shear 


Viscosity  (poise) 


Rate  (sec'1)  DO-1  DO-2  DO-3  DO-4  DO-5  DO-6  DO-7  DO-8  DO-9 


49  0.496 
59  0.504 


0.149  0.495 
0.160  0.489 
0.482 

0.0300  0.159  0.475 
0.0373  0.143  0.478 
0.133  0.473 
0.0300  0.144  0.457 
0.0373  0.139  0.453 
0.0333  0.141  0.460 


0.0366 

0.0359 


0.138 

0.134 


0.422 

0.392 


see  Figures  20  and  21 
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Figure  1:  Pycnometer 
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Figure  4:  Ostwald  Viscometer 
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Figure  5:  Brookfield  Small  Sample  Adaptor 
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Slope  = 0.000199  (g-a/g-cm3) 
Intercept  = 0.9503  (g/cm3) 


(g2/cm3) 


UNCLASSIFIED 


STN  435 


i 


Figure  8:  Dissolved  Density  Determination 
+(See  Equation  2 of  the  Theory) 
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Figures  10  to  18 


Linear  Relationships  Between 
Log  Shear  Rate  and  Log  Shear  Stress 


Log  Shear  Stress  (dynes/cm2) 
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Log  Shear  Stress  (dynes/cm2) 
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Slope  = 0.960 
Intercept  = -0.768 


Log  Shear  Rate  (sec-1) 


+ Non-Newtonian  shear  rates 


Newtonian  shear  rates 


Figure  11:  Dowanol  2 
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Log  Shear  Stress  (dynes/cm2) 
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Log  Shear  Stress  (dynes/cm2) 
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Log  Shear  Stress  (dynes/cm 
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Log  Shear  Stress  (dynes/cm2) 


n 
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Log  Shear  Stress  (dynes/cm2) 
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Figure  17: 


H 1 1 1 1 1 1 1 1 1 F 

1*0  5*0  3.0 

Log  Shear  Rate  (sec-1 ) 


Key 

+ Non-Newtonian  shear  rates 

^ Newtonian  shear  rates 

Dowanol  8 
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Slope  = 0.799 
Intercept  = 1 .346 
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P 

* 


1*0 

e*o  3*o 

Log  Shear  Rate 

(sec-1 ) 

Key 

+ Non- Newtonian  shear  rates 

^ Newtonian  shear  rates 

Figure  18:  Dowanol  9 
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Figure  19:  Effect  of  K-125  on  Viscosity  of  Dowanol 
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Shear  Rate,  \ (sec'1) 


UNCLASSIFIED 


APPENDIX  I 


EXPERIMENTAL  TECHNIQUES 

All  glassware  was  frequently  cleaned  with  a chromic  acid  solution 
and  then  thoroughly  rinsed  with  distilled  water  and  acetone. 

1.  Density  - This  property  was  measured  by  means  of  a pycnometer  with  a 
thermometer. 

The  sample  was  cooled  to  5 to  10°C  below  the  test  temperature 
and  introduced  into  the  clean  pycnometer.  The  thermometer  was  then  inserted 
into  the  pycnometer  taking  care  not  to  trap  any  air  bubbles.  The  pycnometer 
was  then  placed  in  a water  bath  for  10  to  15  minutes  and  allowed  to  come 
to  thermal  equilibrium.  The  excess  fluid  on  the  overflow  arm  was  wiped 
off,  being  careful  not  to  absorb  any  fluid  from  inside  the  arm.  The 
pycnometer  was  then  removed  from  the  bath  and  the  arm  was  capped.  The 
pycnometer  was  wiped  of  excess  water  and  left  standing  to  dry  for  approxi- 
mately 5 minutes,  then  weighed. 

2.  Surface  Tension  - This  property  was  measured  using  glass,  precision- 
bored  capillary  tubing  (see  Fig.  2). 

The  capillary  rise  apparatus  was  set  up  as  in  Figure  2.  The 
test-liquid  was  poured  into  the  test-tube  to  a depth  of  several  centimeters 
in  order  for  the  capillary  tubing  to  be  partially  immersed.  The  test-tube 
was  then  placed  in  a water  bath  and  allowed  to  come  to  temperature  (10 
to  15  minutes).  Pressure  was  then  applied  to  the  liquid  to  raise  the  level 
in  the  capillary,  then  released  to  allow  the  level  to  fall  back  to  its 
equilibrium  level.  Suction  was  then  applied  to  the  liquid  to  lower  the 
level  in  the  capillary,  then  released  to  allow  the  level  to  climb  back  to 
its  equilibrium  level.  Using  a cathotometer,  four  or  five  readings  of  the 
capillary  height  at  equilibrium  were  taken  for  each  test-liquid. 

3.  Vapor  Pressure  - This  property  was  measured  by  the  Ramsey-Young  tech- 
nique using  the  Tobey  modification  (Ref.  1). 
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In  this  apparatus  (Fig.  3)  the  liquid  was  vaporized  from  the 
round  bottom  flask  using  a heating  mantle  controlled  by  an  autotransformer. 
The  vapor  would  condense  on  the  cold  finger  and  flow  down  the  thermometer 
to  a hygrometer  wick  around  the  thermometer  bulb. 

The  temperature  measured  was  that  of  the  liquid  in  equilibrium 
with  the  vapor. 

Equilibrium  was  determined  by  a constant  efflux  time  of  1 drop 
every  5 to  10  seconds  off  the  thermometer  bulb. 

The  pressure  was  controlled  by  a vacuum  pump  attached  to  a 
manometer  and  large  ballast  tank.  The  tank  prevents  sudden  pressure  surges 
from  disturbing  the  equilibrium. 

The  pressure  was  successively  reduced  in  steps,  and  a series  of 
data  points  taken.  The  autotransformer  was  set  in  order  to  maintain  a 
suitable  rate  of  reflux  at  each  reading. 

4.  Viscosity  - This  property  is  unfortunately  more  complex  to  measure 
than  those  previously  described.  Two  different  devices  which  work  on 
different  theoretical  bases  were  used.  The  first  apparatus  is  the  Ostwald 
viscometer  which  works  on  the  principle  of  capillary  flow  (see  Fig.  4). 

The  second  apparatus  is  the  Brookfield  viscometer  which  works  on  the  prin- 
ciple of  coaxial  flow  (see  Figs.  5 and  6). 

a.  Ostwald  Viscosities  - The  apparatus  used  was  the  Canon-Fenske  type 
Ostwald  viscometer. 

To  charge  the  viscometer,  a rubber  tube  was  connected  to  the  large 
arm  and  the  viscometer  was  inverted,  dipping  the  small  arm  into  the 
sample.  Suction  was  applied  to  the  rubber  tube  and  the  liquid  was 
drawn  up  to  the  second  etched  line  to  fill  both  bulbs.  The  visco- 
meter was  then  turned  back  to  its  normal  vertical  position  and  the 
small  arm  was  wiped  clean.  The  viscometer  was  then  placed  in  the 
constant  temperature  bath  and  allowed  to  come  to  thermal  equilibrium 
(15  to  20  minutes).  With  a suction  bulb  attached  to  the  small  arm, 
the  fluid  was  drawn  up  past  the  top  etched  line.  The  suction  bulb 
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was  then  removed  and  the  time  was  recorded  for  the  fluid  to  pass  both 
etched  lines.  Measurements  were  taken  for  times  greater  than  200 
seconds  for  which  kinetic  energy  effects  can  be  neglected  as  recom- 
mended in  the  ASTM  standards  (Ref.  2). 

Three  readings  were  taken  for  each  sample  and  averaged. 

b.  Brookfield  Viscosities  - The  Brookfield  viscometer  was  used  with 
the  Brookfield  Small  Sample  Adaptor.  This  device  uses  the  prin- 
ciple of  coaxial  rotation  with  the  outer  cylinder  stationary. 
Temperature  control  is  achieved  by  a jacket  around  the  outer 
cylinder. 

The  required  amount  of  sample  was  placed  in  the  chamber  and  the 
spindle  was  then  inserted.  The  spindle  and  chamber  were  then  attached 
to  the  viscometer  and  the  sample  was  allowed  to  come  to  temperature 
(10  to  15  minutes)  by  circulating  water  through  the  jacket.  The 
clutch  on  the  viscometer  was  depressed  and  the  motor  turned  on,  rota- 
ting the  spindle  in  the  sample  chamber.  The  clutch  was  then  released 
and  the  dial  was  allowed  to  rotate  until  the  pointer  stabilized  at  a 
fixed  position  on  the  dial.  The  lowest  possible  speed  to  obtain  a 
reading  was  determined  then  further  readings  were  recorded  from  the 
lowest  to  the  highest  possible  speed.  Torque  was  then  evaluated  from 
these  readings. 

Two  slight  modifications  were  built  into  this  device.  A longer 
spindle  and  chamber  were  constructed  to  determine  if  end  effects 
were  present  in  the  smaller  apparatus  and  to  study  liquids  with  lower 
viscosities.  Also,  a thicker  spindle  was  constructed  for  the  small 
device  in  order  to  extend  the  range  of  shear  rates. 
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APPENDIX  II 


RHEOLOGY  REVIEW 

Rheology  is  the  study  of  the  reaction  of  materials  to  applied 
forces.  This  brief  review  is  based  extensively  on  the  material  in 
Middleman's  book:  The  Flow  of  High  Polymers  (Ref.  4). 

There  are  three  basic  categories  used  to  describe  materials. 

These  are: 

i)  Hookean  Elastic  (ideal  solid) 

ii)  Newtonian  Fluid  (ideal  liquid) 

iii)  Anything  Else 

The  first  two  categories  are  simple  mathematical  models  which 
seem  to  describe  many  materials  under  normal  atmospheric  conditions.  The 
third  category  consists  of  a multitude  of  mathematical  models  each  claiming 
to  represent  a particular  class  of  materials  which  do  not  fall  into  either 
of  the  first  two  categories.  These  rheological  models  (mathematical 
equations)  are  useful  to  help  predict  how  a material  will  react  under  new 
or  different  environmental  conditions. 

The  rheological  equations  take  the  form  of  relating  the  stress 
tensor  (internal  forces  acting  over  an  area)  to  the  strain  tensor  (the 
actual  displacement  of  the  material  relative  to  an  initial,  relaxed  state). 
Some  of  the  more  commonly-used  shear  relations  include: 

i)  Hookean  Solid  (elastic): 

Tij  = Gi 

stress  « strain 

ii)  Newtonian  Liquid  (viscous): 

Tij  = V 

stress  « strain  rate 
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iii)  Power-Law  Fluid: 

T. . = n Yn 
ij  oT 

iv)  Bingham  Plastic: 


T1j  = < ’o> 


1J 


To  + V-"-(tU  > To> 


v)  Viscoelastic  Fluid  (Maxwell  Model,  Linear): 


. no  6tij  . . 

x . . + — = n y 

1J  G St  0 


vi)  Viscoelastic  Fluid  (Kelvin  Model,  Linear): 

Tij  = V + 6y 

vii)  Viscoelastic  Fluid  (Spriggs  Model,  non-Linear): 


, o \ _ 

+ — (F  T.,)  - n Y 


1J 


1J 


. ...  where  F is  a special  derivative. 


The  above  seven  relationships  represent  attempts  to  model  real  phenomena 
by  equations  of  various  degrees  of  complexity.  Many  other  models  exist. 
For  some  models,  the  constants  can  vary  with  time  in  order  to  simulate 
thixotropic  or  rheopectic  fluids.  Modern  models  use  tensor  notation 
thereby  making  the  equation  independent  of  any  coordinate  system. 

When  trying  to  fit  experimental  results  to  a rheological  model, 
the  flow  situation  must  be  kept  simple.  There  are  two  basic  types  of 
experiments  used  to  determine  rheological  coefficients  for  fluids: 


1 . Simple  Shear  Experiments 

The  coefficient,  n,  is  determined  by  dividing  the  shearing 
stress,  T^j,  by  the  shearing  strain  rate,  y.  n may  be  a function  of  y. 
Three  of  the  most  commonly  used  techniques  are  illustrated  in  the  following 
sketches. 
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i)  Coaxial  (Couette)  Flow  (4) 


Y = 2q  I s2'3mm 
p=o 

Tr  = T/Utt^L) 

where  m evaluated  at:  t = s2|\ 


ii)  Capillary  (Poiseuille)  Flow  (4) 
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2.  Elongation  Experiments 


The  coefficient,  ng,  is  determined  by  dividing  the  normal  stress, 
t .j  • » by  the  elongational  strain  rate,  y . Very  few  experiments  can  be 
performed  with  low  viscosity  liquids  in  this  area.  Some  examples  are: 

i)  Simple  Stretch  (4) 


U 


Ye  - U/L 

F/  ( ttR2  ) 


! 


Viscous  Liquid 
- Clamped 
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Complications  Due  To  Viscoelasticity 

Many  liquid  solutions  consist  of  long  polymer  chains  which  are 
incorporated  to  give  a proper  "texture"  to  a fluid.  These  chains  act 
somewhat  like  tiny  springs  which  can  be  stretched  by  a straining  force 
and  will  recoil  when  that  force  is  removed.  These  molecules  are  res- 
ponsible for  the  viscoelastic  properties  of  the  fluid. 

Viscoelasticity  is  commonly  modelled  by  the  two  types  of  linear 
equations  described  earlier.  These  linear  equations  assume  that  the 
elastic  and  viscous  components  can  be  treated  separately.  They  are  suit- 
able for  modelling  relaxation  phenomena  but  they  cannot  predict  normal 
stresses  which  occur  during  shearing  flows.  In  this  situation,  a non- 
linear equation  (such  as  the  Spriggs  model)  must  be  employed. 

Viscoelastic  models  require  the  determination  not  only  of  the 
viscosity  coefficient  for  flow  but  also  the  elasticity  coefficient  for 
elastic  stretch  and  recoil.  Considering  the  Maxwell  model  (commonly  used 
for  low  to  medium  viscosity  liquids),  for  an  abrupt  change  in  strain  rate 
the  stress  relaxes  according  to: 

Tjj  = -Qexp  (-Gt/n) 

Therefore,  a plot  of  In  t. . vs  t would  yield  G/n  and  with  n determined  from 

^ J 

a simple  shear  experiment,  G could  be  calculated.  Another  technique  is 
to  apply  a sinusoidal  strain  rate.  A dynamic  viscosity  coefficient  is 
then  defined  by: 

n*  = n'  - in"  = TQsin6/(wYo)  - i t0cos6/(wy0) 

...  where  6 is  the  phase  angle  between  stress  and  strain.  The  dynamic, 
elastic  coefficient  is  then  given: 

G*  = iw  n* 

For  the  non-linear  viscoelastic  models,  one  needs  to  measure  normal  stresses 
which  occur  due  to  shear  flow.  This  defines  a "cross-viscosity"  phenomenon. 
Fortunately,  there  are  instruments  which  are  capable  of  doing  this,  such  as 
the  famous  Weissenberg  cone  and  plate  rheogoniometer. 
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F,  W = forces  (F) 

G = Young's  elastic  modulus  ( F/L2 ) 

L = length  (L) 

m = slope  of  log  fi  vs  log  t 

aP  = pressure  drop  (F/L2) 

Q = volumetric  flow  rate  ( L 3 1~ 1 ) 

R,  r = radius  (L) 

t = time  (t) 

T = torque  (FL) 

u,  v = velocities  (L/t) 

w = frequency  (t-1) 

x = distance  (L) 

i|>  = angle 

6 = phase  angle 

nQ  = coefficient  of  shear  viscosity  (Ft/L2) 

ng  = coefficient  of  extensional  viscosity  (Ft/L2) 

y = shear  strain 

y = shear  strain  rate  (t_1) 

o = rotational  speed  (t-1) 

t = shear  stress  (F/L2) 

= normal  stress  (F/L2) 


r ’ 
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